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This  paper  presents  a  nev  model  for  the  machine  calculation  of 
fallout  patterns,  ihe  model  is  divided  into  three  parts!  (1)  meteorologl*- 
cal,  vhirh  rnnsiderfi  Hind  velocity  and  the  fall  vt^lucity  of  the  particles, 
(2)  geometric,  vfilch  concidern  distrlhution  of  activity  on  paxticlGS  in 
the  cloud  as  a  function  of  height  and  particle  size/  and  (3)  radiological, 
which  considers  the  radiation  doses  of  particles  on  the  ground.  The  use 
of  the  model  is  illustrated,  and  the  changes  in  the  model  that  will  pro¬ 
vide  an  optimum  model  (currently  hoing  develoned)  are  d^J  RcuAfind . 
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A  now  model  for  the  nvichlno  calculation  of  faLlofut  dencrlbed  in  this  Research 
Memorandum  la  the  latest  reoult  of  reoGarch  which  began  at  RAND  In  1952.  An 
earlior  model  to  explain  the  obocrved  facta  of  radioactive  fallout  vaa  described 
in  other  publicationo  (liM-l855j  RM-1371#  R-265-AECj  and  )•  The  present 

Research  Memorandum  deocribea  a  new  and  more  flexible  computational  model  to 
Inveatlgate  the  effect  of  variations  of  input  data  and  to  discover  the  effect 
of  parameters*  The  model  io  divided  Into  three  parts:  meteorological,  geo¬ 
metric,  and  radiological. 

The  initial  test  of  the  model  was  made  with  the  fra^^ntory  Castle-Bravo 
wind  data.  The  feaoihility  of  the  approach  has  been  shown,  and  the  accuracy 
of  the  model  within  the  limits  of  the  available  data  has  been  verified.  The 
next  set  of  data  processed  can  proceed  on  the  aseuiuptions  that: 

1.  small  changes  in  the  distribution  functions  will  produce 
only  negligible  changes  in  the  pattern; 

2.  the  wind  patterns  ore  sufficiently  accurate  to  put  limits 
on  the  possible  distribution  functions; 

3.  the  gross  effects  of  the  fallcmt  from  large-yield  shots 
are  not  affected  by  stem  material. 

It  may  be  Inferred  that  moat  of  tho  radioactive  material  la  lodged  on  particles 
in  a  very  narrow  range  of  sizes  and  that  it  Is  concentrated  into  a  very  narrow 
range  of  heights  at  the  time  of  stabilization. 

NOTE;  This  Research  Memorandum  is  classified  SECRET-RESTRICTED  DATA.  This 

In  Brief  is  classified  SECRET  and  does  not  contain  any  RESTRICTED  DATA. 
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I,  CTfUODUCTION 

The  procecs  of  dcvelopirv-^  a  model  Dmt  vlll  adeqiiatcly  explain  the 

observed  facte  of  radioactive  fallout  van  begun  at  HMD  in  the  fall  of  1952. 

Thu  pii^Biuul  ile^cribod  In  reference  (1)  vao  the  rceult  of  thoce 

early  inveatigatioua*  The  sin^ile  mathoTnaticol  modcl^^^  and  the  Bubeequent 

(2) 

work  of  other  invest Igatoro'  cerved  to  confirm  the  validity  of  the  approach 
to  the  problem,  Proccdureri  were  devised,  for  hand  computation  and  machine 
coii5)utatlon  which  further  confirmed  the  bade  principles  but  also  served 
to  shov  that  there  were  deficiencies  in  the  input  data.  To  Investigate 
adequately  the  effect  of  variations  of  input  data  and  to  discover  the  effect 
of  variation  of  parameters  a  nev^  more  flexible  coc^iutational  model  was 
needed.  This  report  deals  with  the  development  of  such  a  model. 

Before  presenting  the  development  of  the  new  model  it  may  he  profitable 
to  review  briefly  the  physical  concepts  and  the  first  machine  model.  To 
describe  the  phenomenon  of  falloutj  we  must  uegin  with  the  nucleat-  detoua- 
tion.  TiiC  diGtinguiching  feature  of  s  detonation  of  this  type  is  a 

tremendous  release  of  cnerf*y  almost  Instontaneoualy  and  within  a  very  small 
space.  This  release  of  cner/^y  creates  a  hot  bubble  of  air  which  rises  be¬ 
cause  of  its  bouyoncy.  The  characteristic  muGhroQm-*shapcd  cloud  associated 
with  nuclear  explosions  is  formed  by  water  vapor  condensed  as  a  result  of 
vertical  lifting.  The  temperature  grfidlcnt  between  the  hot  bubble  and  the 
surrounding  atmosphere  causes  a  strong  toroidal  circulation. 

Fission  fragments,  together  with  any  other  debris  mixed  into  the  fire¬ 
ball  at  an  early  stage,  are  carried  aloft,  and  as  the  cloud  rises  Bome  of 
the  material  may  be  left  behind,  forming'  a  wake  or  stem,  Tlie  cloud  will 
stop  rising  when  the  energy  available  to  lift  it  has  dissipated.  At  this  time 
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(from  4  to  6  minutes),  which  le  referred  to  aa  the  time  uf  ^tablll^ation, 
the  apatlal  diatrihutlon  of  the  radioactivity  provides  one  aet  of  the  Initial 
condltlona  required  In  the  fallout  model# 

Within  this  stabilized  cloud  are  virtually  all  the  fiesion  products 
formed  by  the  nuclear  detonation  and,  in  some  coses,  by  induce^d  activity. 
Radioactive  particles  may  be  formed  by  condensation  of  vaporized  materleil 
or  by  the  coalescence  of  fission  products  with  solid  or  liquid  earth  parti-*- 
dee.  The  size  of  the  paxticlee  is  determined  largely  by  the  initial  con¬ 
ditions  of  the  burst;  if  little  solid  material  Is  token  Into  the  cloud,  the 


vlll  in  lihcllhccd  be  crdl;  if. 


— kgr,  VjLctvacC  xt> 


detonated  on  or  near  the  surface,  many  tons  of  earth  will  be  swept  into  the 
cloud,  and  a  sizeable  fraction  of  the  fission  products  will  be  contained  in 
large  particles.  Tlxe  size  of  the  particles  containing  radioactive  elements 
has  an  Important  bearing  on  the  fallout  problem.  The  radioactivity  located 
In  or  on  large  particles  will  be  brought  rapidly  to  the  ground,  while  that 
lodged  In  small  particles  will  remain  suspended  in  the  atmosphere  for  long 


periods  of  times  Another  of  the  initial  conditions  required  In  the  fallout 
model,  therefore,  la  the  distribution  of  activity  with  partlclo  size  within 
the  stabilized  cloud* 

Given  the  stabilized  cloud  and  the  associated  distribution  of  activity 
with  particle  size  and  in  space,  it  la  necessary  to  determine  the  position 
of  the  active  particles  after  they  fall  to  the  ground.  The  rate  of  feQl  is 
of  primary  luportance  here,  and  on  aerodynamic  law  of  fall  which  takes  Into 
account,  by  raeeins  of  drag  ceeff Iciants,  inertial  as  well  aa  viscous  forces 
was  used. 

As  the  particles  start  to  fall  through  the  atmosphere,  they  will  be 
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carried  by  the  vind.  There  are  a  great  many  wind  effects  which  have  to  be 
taken  into  account.  First,  the  whole  cloud  will  be  moved  bodily  by  the  mean 
wind  at  the  cloud  level.  At  the  same  time,  there  will  be  shearing  effects 
due  to  the  variation  of  the  wind  with  altitude,  and  there  will  be  turbulent 
effects,  small  scale  eddleo,  which  will  tend  to  Bpread  and  tear  the  cloud. 
The  neceeaory  elements  of  a  fallout  model  may  be  summarized: 

1.  Distribution  of  activity  with  space 

2.  Distribution  of  activity  with  particle  size 

3.  Fall  velocity 
Wind  structure 
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II «  FIICT  MACmWK  mi)V2j 

The  first  machine  model  need  at  RAND  is  reported  in  detail  In  reference 
(4).  The  cloud  was  ascumed  to  be  a  cylinder  with  a  uniform  dietributlon  of 
activity  In  the  horizontal  (Fig.  1)  and  on  exponential  decrease  frcHn  the 
base  of  the  cloud  upward  (Fig*  P.)*  This  decrease  with  height  vaa  predicated 
on  the  aaBU^^)tlon  that  the  number  of  active  partlclco  per  unit  ma^  of  air 
is  constant  in  the  munhrocm.  The  ixniforra  horizontEil  distribution  is  the 
most  Bljr5)le  of  all  poeelble  dlstributlonG .  The  distribution  of  activity 
with  particle  size  is  baaed  on  a  dlDtribution  which  reproduced  the  fallout 
pattern  of  the  Jangle  Surface  shot  (Fig*  3)* 

In  the  model,  port  idea  were  chosen  to  represent  the  fallout  at  every 
1000  ft  of  elevation  at  the  center  of  the  cloud  for  values  of  particle  size 
which  divided  the  activity  size  distribution  into  100  equal  Intervals.  Thus 
for  each  elevation,  each  particle  size  represented  the  same  amount  of  actlv^ 
ity.  The  turbulent  effects  were  ignored  and  the  ground  position  of  each 
particle  from  each  elevation  was  assumed  to  represent  the  center  of  a  circle 
whose  radius  was  equal  to  the  radius  of  the  stabilized  cloiv^e  The  activity 
at  a  point  vaa  calculated  by  adding  the  contributions  from  those  circles  which 
contained  the  point.  This  syatem  of  ca3*culatlon  is  fairly  efficient  insofar 
as  machine  time  is  concerned;  however,  the  method  requires  a  complete  recoding 
for  a  change  in  the  model  because  the  meteorology  is  connected  to  the  geomet¬ 
ric  and  particle -*Bize  parameters  through  the  choice  of  particle  sizes  to  be 
used.  The  new  coo^jutation  scheme  should  permit  variations  of  scane  parts  of 
the  model  without  changing  others. 

Another  feature  of  the  first  ccm^iutation  scheme  was  the  Incorporation 
of  the  decay  scheme  into  the  procedure.  The  activity  was  assumed  to  decay 
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Dimensions  of  cloud 


Width  (n  mi) 

Fig.  I  — Assumed  geometry  of  a  cloud  resulting 
from  a  large-yield  nuclear  detonation 
(As  used  in  the  first  RAND  machine  model) 
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Fig. 2 — Assumed  distribution  of 
radioactivity  in  space  following  a 
large-yield  nuclear  ueiunuiiun 

(Ag  used  in  the  first  RAND  machine  model) 
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3 — Cumulative  distribution  of  radioactivity  with 
particle  size 

(As  used  in  the  first  RAND  mactiine  model; 
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~1  2 

vlth  the  t  *  laVj  and  the  result  a  vere  tranafonced  to  doae  rates  and  Inte^ 
grated  doaca.  While  this  achomo  worked  well  for  fission  devices,  the  vari¬ 
ability  introduced  by  clean  and  salted  weapons  makea  it  li^portant  to  separate 
the  radiological  effects  from  the  calculations  of  particle  deposition. 
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III,  TIIKWKW  MODI-JL 


A,  General 

In  ordor  to  Incnre  flexibility  In  the  new  method  of  confuting  fallout, 
the  problem  was  separated  into  three  con^Donent  parts.  The  first  port,  called 
the  meteorological  port,  is  the  determination  of  the  positions  on  the  ground 
of  a  fixed  set  of  particles  atorting  from  fixed  heights  In  the  atmosphere. 
There  ore  only  two  sots  of  the  basic  variables  introduced  In  thlo  section 
of  the  commutations,  1.^  ,  the  fall  speed  and  the  velocity  of  the  wind. 

The  second  section,  which  is  referred  to  as  the  geometric  part,  is  a  method 

for  i-nM  Kill  i  mi'm  mh  m  I'tUiCvioU  OI  XHC  0.0*10 xuy  Liie 

bilized  cloud,  using  the  ground  distribution  of  particles.  This  section 
involves  the  distribution  of  activity  with  height,  radius,  and  particle 
size.  The  third  part,  which  is  called  the  radiological  part,  deals  with  the 
ELnalysls  of  the  radiation  dosof^es  that  ore  received  from  the  distribution 
of  activity  on  the  ground.  Due  to  the  variable  nature  of  the  radiological 
products  no  attenmt  will  be  made,  in  this  model,  to  confute  radiological 
effects  on  the  ground. 


B.  The  Meteor ologlc  Calculations 

The  ground  position  of  a  particle  falling  with  a  velocity  w(r,h)  through 
a  wind  field  \V(x,y,h,t)  is  given  by 


(D 


v(r,h} 


* 

In  contrast  to  the  first  model  which  was  run  on  the  n3M“702,  the  second 
model  was  programmed  for  and  run  on  the  IBM-Y04  with  on  augmented  memory. 
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For  numerical  work  thia  can  be  expreafied  ae  a  Bura  and  divided  into  tvo  codT' 
ponent  parte 


X 


m 

.  V 

ir-O 


v^(r>h) 


Ah 


1 


y 


m 

\ 

/ 

i:.*n 


(r/hl 


Ah, 


The  fall  vclocitiea  ar'o  rnmputed  by  the  aerodynamic  methods  outlined 
in  the  references  mentioned, A  detailed  onalyslc  of  the  problem  of 
rail  velocities  will  be  the  oubject  of  another  report,  It  is  aufiiclent 
here  to  state  that  aerodynamic  x’ates  for  spheren  have  be^n  ctsed  In  accordance 
vith  the  results  of  the  experimenta  reported  in  reference  (5)*  The  values 
of  the  fall  speeds  ore  presented  in  Fig.  4, 

It  has  been  shown^^^^^^  that  the  time  and  space  variations  of  the  vind 
are  Important  in  predicting  fallout^  and  therefore  the  model  must  be  able  to 
account  for  time  and  space  variations#  The  wind  field  Is  really  a  three- 
dlmenBloned  vector  field  which  Is  continuous  in  time  and  epace#  To  properly 
take  this  into  account,  continuous  functions  should  be  used#  Uhfortunately 
measui'ements  of  the  wind  field  give  discrete  values  of  the  horizontal  conpo- 
nent  at  discrete  times.  The  vertical  coirponent  cannot  be  directly  measured 
by  any  instrument,  but  it  is  known  to  be  small  ccaipared  to  the  horizontal 
cnnmor  jiiu,  lac  vux  uumxxouunt  will  be  neglected  because  it  cannot,  vith 

current  techniques,  be  calculated  vith  sufficient  accuracy. 

The  eaat-^wcBt  and  nortlr-couth  componentB  of  the  horizontal  wind  can  be 
measured  and  vlU  be  used  separately*  Thus,  if  represents  one  of  these 
conponents.  It  con  be  moacui^cd  at  discrete  intervals  of  space  and  time#  In 
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order  to  CBtimate  the  vlnd  nffectinti  the  particle  at  any  point  and  at  any 
instanti  a  four  dimensional  linear  intci’polation  ocheme  vlU  be  used. 

Suppose  at  time  t,  the  particle  is  at  a  height  H,  at  the  position  IJ. 
Let  represent  the  velocity  ccii5)onont  in  the  a  direction  at  the  point 

ij,  the  height  k,  and  the  time  n.  If  there  ore  three  points  vhere  the  wind 
1b  known  at  two  timea  and  at  intervals  of  5000  ft,  a  first  interpolation  can 
be  laode  In  ele'^'stlon: 

Hq  - 

^orijon  ^odn-ln  ^  X 


The  next  interpolation  is  for  position  and  is  given  by 

,  r 

^txoo  "  ^a33  *  D  ^  ^  ^^1  -  ^2^  ('^«33  ”  ''«11^  -  ^^1  “  ^3^  ^^a22  “ 

+  y  (V^2  ”  “  Xg) 

+  (Yi  -  y^)  (v^  -  v^)  -  (y^  -  y^)  -  V^) 

T  y^  (x^  -  xg)  (V^33  -  -  (x^  -  X3)  (V^2  " 

where  D  =  (x^  -  x^)  -  y^)  -  (x^  -  x,)  {y^  —  y^)  and  "...on”  is  assinned 

to  be  added  to  the  subscripts  of  all  the  V*e. 

Tlie  fimil  interpolation  is  in  time  and  is  given  by 


^nOOOO 


V  + 

nOOO-1 


(V 


nOOO-1 


—  V  ) 
QOOOl^ 


(3) 


If  this  is  done  for  both  con^onents  x  and  y,  then  the  velocity  at  the  posi¬ 
tion  of  the  particle  is  determined. 

With  this  interpolation  process^  a  particle  is  stai'ted  from  fO'ound 
zero,  the  vlnd  componento  arc  computed  from  Kqs.  (l),  (2),  and  (3).  These 
winds  ore  multiplied  by  a  thicknosG  Ah  and  divided  by  the  mean  fall  velocity 
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In  the  layer;  thie  ylcldn  a  new  ponitlon  (x^y).  The  above  process  is  re¬ 
peated  with  this  new  position  and  a  nt.‘W  chon'^e  In  position  is  con^juted. 

The  process  is  repaated  until  the  particle  reaches  the  pyround*  A  total,  of 
7500  particle-height  ccsnbinatlonc  are  traced  to  the  ground  in  this  manner. 

The  time  required  for  a  r.lven  pai’tlcle  to  reach  the  surface  from  a 
given  height  is  also  computed  from  the  fozTnula 

d  “  /  w(r,hy'  (!)•) 

The  ground  position  and  tlme-^ovn  for  each  particle  from  each  height 
is  stored  on  a  tape  so  that  the  wind  system  con  be  used  to  cheek  variations 
in  geometry  which  might  arise  from  variations  in  yield  or  burst  conditions* 

C*  Errors  in  Meteorological  Port 

It  is  perhaps  worthwhile  to  form  some  sort  of  estimate  of  the  accuracy 
which  may  be  obtained  in  the  meteorological  part  of  the  problem.  In  order 
to  make  this  estimate  a  great  many  sln^ilif ications  will  be  mode  so  that  the 
estimates  will  yield  only  the  order  of  magnitude  of  the  wind  error. 

There  ore  several  major  asoun^tiono  which  need  to  be  made  in  this 
analysis: 

1.  Wind  errors  6(o;CT^)  arc  distributed  according  to  a  circular 

2 

normal  distribution  with  zero  mean  and  a  variance  of  • 

2.  The  variance  of  the  errors  io  not  a  function  of  elevation. 

5.  The  correlation  of  errors  is  a  function  only  of  the  separation 
between  winds  and  decreases  rapidly  with  height  difference. 

Us  For  the  purposes  of  this  error  analysis  only^  the  fall  velo¬ 
cities  will  be  approximated  by  W  -  W^(r)e^. 
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The  fall  velocitioo  are  correct,  and  no  vertlcaLl  air  motion 


Nov  the  calculated  position  of  a  given  particle  vlll  be  its  actual  po¬ 


sition  plus  the  Integral  of  the  error  term 


h  e  “^4  +  y  ■  wfrVhT 


/''cording  to  assun^ition  1,  =  o,  and  E  «  o  but  a_  must  be  calculated  from 

£i 


H  H  fz — 

2  C  C  '  1  2  ^ 

""  “  \  ,1  'wThTTwI 


dh.  dlu 

J.  d. 


with  asBun^jtiou  2  above  this  becoiaee 


2  2  /*  (^ 

‘'e  “  “e  j  ^  VTh^TwTi^r 


where  correlation  of  errors  with  height.  Nov  in  accordemce  with 

aB3un5)tlon  5/  ^  where  p  la  of  the  order  of  l/2.  This  means  that 

the  correlation  ia  less  than  0.1  for  separations  greater  than  5000  feet. 
According  to  assun^itlon  4,  W(r,h)  =  W^(r)  where  the  largest  vslIucs  of 

a,  for  500p  particles,  lo  about  0.016.  Thus  p  in  an  order  of  magnitude 
greater  than  nr. 

Inserting  these  values  in  Eq.  (7) 


■  -^rT 

rt  o'  ' 


-P(lL  -  hg  ) 


where  the  subscript  r  has  been  added  to  the  notation  for  the  error  variance 
to  indicate  that  r  is  treated  as  a  parameter.  By  making  a  change  of  variables 
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Eq.  (8)  bccomoB 


?  E 

"rE  =  -OTT 


n  ^  4 


-t 


‘  (pt  4  ra)  , 

'  do  dt 


(9) 


Which,  \^on  inte(;i'rition,  bcccmuj; 


<(r)  ) 

"  I 


(l-e 


(p-a)ll^ 


J 


Now  if  p  is  of  the  order  of  OO  find  nr  la  of  the  order  of  0.01,  and  If  H  Is 
greater  than  10  kilofect,  the  temia  e  '  '  ^  and  e  '  ore  small.  con5)ared 

to  1  and  p-CT  and  pio  are  very  neoi-ly  equal  to  p.  Therefore,  for  H  >  10 
kllofoet 


4  s  —4--  (1  - 

op  (r) 
o 


(11) 


Equation  (11)  vao  used  to  evaluate  the  ratio  of  to  for  a  few 
values  of  r  and  II  in  ordf'r  to  provide  a  means  of  estimating  the  expected 
errors  in  a  fallout  pattern.  'I'liosc  oi'c  shown  in  Table  1. 


Table  1 


Ratio  of  o 


I'K/ 


ns  a  EuiietJ-on  of  llcirlit  and  Particle  Radius 


H 

1* 

(|l) 

Xkf) 

2h 

(Jo 

151 

600 

10 

3.101 

.7't7 

.2)^3 

.003 

20 

ii,l6o 

1.02.3 

.326 

.OBk 

30 

5.000 

1.190 

.375 

.095 

*10 

5.633 

1.321 

.iiir? 

.103 

6o 

6.573 

1.306 

.I1&’ 

.11*1 

00 

7.  OTP 

1.6,28 

.I191 

.118 

100 

7.303 

1.701 

.51? 

.118 
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D«  Gi^omctrlc  Vart 

The  meteoroloKlc  part  of  tlic  problem  provides  a  method  for  determining 
the  ground  position  of  m/iny  poi'ticlt'u  originating  on  many  helfdrts  over 
ground  zero.  If  an  nihlitionni  Mnuitn^^tlon  Is  made,  namely  seme  lav  for  the 
distribution  of  poi-ticlea  in  the  horizontal,  the  ground  positions  of 
many  particles  from  many  positions  in  the  cloud  could  be  determined.  If 
a  fraction  of  activity  per  unit  volume  per  unit  particle  size  can  be 
assigned  as  a  function  of  position  in  the  cloud,  then  it  is  possible  to 
deterxo^ne  the  fraction  per  unit  area  in  the  fallout  pattern. 

The  first  slu^jliflcation  that  vlll  be  made  is  to  assume  that  the 
distribution  is  radijilly  symmetric  about  the  ground  zero  pouitiun.  A 
second  simplification  is  to  assume  that  tm'bulent  diffusion  is  not 
iuroortont,  so  tliat  a  pai'tlcle  at  a  distance  R  from  the  ground  zero  in  the 
cloud  would  strike  the  ground  at  a  distance  R  from  the  position  of  a 
particle  of  the  same  size  oi iginaling  over  ground  zero  in  the  fallout 
pattern. 

Let  A(R,H,r)  be  the  fraction  of  the  device  per  unit  voltnne  per  unit 
particle  size  In  the  cloud  at  stabilization.  Since  A  is  a  fraction  or 


proportion 

zTt  R  H 

max  max 

r  r 

r 

max 

r 

J 

5  5  1 

A(R,H,r)  dR  dll  dr  de 

(12) 

000 

0 

Since  the  function  of  three  vai'inblec  is  rat)ier  Intractable  it  will  also  be 
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aosumed  that 

A(R,H,r)  -  Aj^(r)  A^(H) 

To  dotormine  the  surface  density,  conalder  an  arbitrary  point  y^) 

in  this  ffillont  pattern.  Near  this  point  P  there  will  be  found  f.bp  position 
of  a  particle  8l7.e  r^^,  lallln//  from  a  height  call  this  point  y^) 


With  the  assunqption  that  the  cloud  conQ)oaed  of  particles  of  size  r^ 
is  not  distorted  by  turbulence,  the  particles  of  size  r^  from  height 
that  land  at  P  will  come  from  a  point  ■=y  ^  (y^-yo)^  distant 

from  the  center  of  the  cloud.  The  volume  concentration  due  to  this 
particle  size  and  this  height  at  the  point  P  is  therefore  given  by 
Ai(ri)  A^CR^)  A^(R^).  Nov  in  a  similar  vay  the  pai’tlcle  size  r^  may  fall 
from  a  height  at  the  point  ^  volume  concentration  of 

Aj^(r^)  AgCHg)  A^(Ry) 
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Nov  all  of  the  particles  from  to  Hp  will  land  at  P  with  radii  lying 
between  find  do  tfint  the  lute.-Tation  of  t)ie  activities  as  determined  by 
the  prodiictc  A^ApA^  between  11^  anil  R,^  will  give  the  nurfnce  density  of  Rctiv^ 
ity  at  P  from  partieleo  of  nli'e  r^.  A  oimilor  argument  for  another  particle 
Bl7.o  will  yield  n  ourfeco  density  duo  to  partieleo  of  oize  If  the 

integration  io  performed  over  particle  size  from  to  r^_,  the  result  is 

JL  t: 

the  fraction  of  the  device  per  unit  area  contributed  by  all  particles 
tween  and  r^  from  heights  between  and  Hp.  If  the  integration  is  ex¬ 
tended  over  all  r  and  all  H,  the  total  fraction  down  at  P  la  obtained.  The 
fraction— down  per  square  kilo loot  Is  calculated  lor  a  set  oi  pre- -chosen  grid 
points  P. 

It  is  IntcreDtlng  to  note  tbnt  the  geometry  of  the  cloud  is  expressed 
by  the  two  distribution  functions  A^(R)  and  A^(}l).  At  large  distances,  A^ 
reduces  to  zero  and  no  contribution  is  made;  this  continuous  distribution 
is  contrasted  with  the  cooky^cuttcr  technique  of  the  first  RAND  model.  The 
errors  Introduced  thi-ouvdi  the  distribution  functions  can  best  be  handled  by 
a  variation  of  partimctcr  techniques.  Thus  several  runs  may  be  made  using 
different  "A”  functions  with  the  some  wind  patterns,  in  order  to  estimate 
the  accui’acy  of  the  gecpnctrlc  part, 

E,  Preliminary  Estljnatco  of  the  Distribution  Functions 

In  on  attempt  to  delineate  the  distribution  of  the  radioactivity  In  a 
stabllzied  cloud,  rockets  were  eliot  through  some  of  the  Redwing  clouds  at 
7  minutes  after  deiicnation.  The  initial  estimates  of  kiloroentgeno  per 
hour  at  7  minutes  was  published  in  ITR-'135^*^^^  It  is  hoped  that,  when 
final  reports  are  prepared,  all  of  the  data  will  be  useful.  Until  such 
time,  the  Zuni  data  will  be  ur:.cd  to  obtain  a  preliminary  estimate  of  the 
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dlotrlbutionn . 

Tlio  '/ami  vookt?Ly  ptiuiunl  very  nf'arly  ovt*r  Kcound  7.(?ro  and  therefore  vere 
aHBinned  to  reprooent  n  vertical  plane  throw-di  'the  dianieter  of  the  cloud.  In 
order  to  further  olmiillfy  the  problem  circular  oyrametry  about  pr^^ound  zero 
will  be  asoTuned.  The  data  were  then  plotted  at  the  appropriate  height  and 
dlDtance  from  Jp^ound  zero  and  nmooth  contours  of  the  reported  rocket  read¬ 
ings  vere  drawn  (t'lg.  5)* 

From  the  smoothed  contoui*u,  the  total  activity  In  the  cloud  was  coc^^uted 
by  a  numerical  integration.  From  the  total  activity  and  the  smooth  contours^ 
the  fraction  of  activity  per  cuuic  Klloi'oot  was  cui^uted  aa  a  fimeticn  of 
height  and  radial  distance  is  Dhf)wn  in  Table  2. 

The  estimate  of  A^(H)  was  clioscn  as  the  marginal  distribution  which  is 
obtained  from  the  numerical  integration  of  R  A(R,H)  dR  at  several  different 
heights.  Tills  distribution  is  shown  in  Fig.  (6).  The  estimate  of  A^(r) 
was  chosen  as  the  integration  of  A(R,H)  dH  for  several  different  values  of 
R.  Figure  (7)  shows  this  distribution. 

The  distribution  of  activity  with  particle  size  was  chosen  os  the 
distribution  presented  in  reference  (1).  Accumulative  distribution  was 
plotted;  points  vere  rood  from  the  cumulative  function;  and  numerlceG.  differ^ 
entiation  was  used  to  determine  values  of  A^(r).  Since  this  lb  a  badly 
skewed  distribution,  it  la  convenient  to  use  as  a  variable  ln(r)  instead  of  r. 
Figure  (8)  shows  the  activity  as  a  function  of  size  for  this  distribution. 
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Fig. 6  —  Smoothed  distribution  of  activity  with 
height  from  Zuni  rocket  data 
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Tlie  Initial  teat  of  tlio  vliich  vill  be  referred  to  ac  "run  Ij" 

vaa  made  vith  the  Cantle— Bravo  wind  data  taken  from  Dean  and  Ohmatead.  The 
detail^a  are  presented  in  the  Appendix. 

The  ground  position  of  a  few  selected  particles,  which  were  read  from 
the  results  of  the  meteorological  port  of  the  confutation,  are  shown  in 
Fig.  {9)f  which  also  shows  the  circular  standard  errors,  derived  from  the 
discussion  of  meteorological  errors,  vith  the  assunftlon  of  on  error  stan- 
dord  deviation,  o^  ,  equal  to  5  knots.  It  appears,  from  this  figure,  that 
random  errors  of  this  magnitude  will  not  produce  large  errors  in  the  orien¬ 
tation  and  general  configuration  of  the  pattern.  It  does  seem  reasonable 
to  assume,  however,  that  this  type  of  error  may  cause  a  considerable  error 
in  the  fraction  of  fallout  at  some  point  near  the  boundary  of  the  fallout 
area.  It  must  be  enfhaoized  that  these  wind  error  estimates  do  not  include 
any  gross  error  in  definini^  tlie  vector  wind  field. 

Several  variations  of  the  activity  functions  were  tried  with  the  basic 
wind  pattern  of  run  1.  The  purpose  of  these  calculations  was  to  make  some 
qualitative  estimates  of  the  effect  of  changing  the  distribution  functions. 
The  details  of  these  calculations  are  given  in  the  Appendix  and  are  referred 
to  as  run  la,  lb,  etc.  The  inferences  which  may  be  drawn  from  an  inspection 
and  coo5)arison  of  the  results  of  the  calcvilationa  will  be  taken  iq)  in  order. 

The  difference  between  runs  la  and  lb  was  only  in  the  function. 

The  function  used  in  lb  wo.s  a  lo/^normol  distribution  with  the  same  mean 
and  standard  deviation  as  the  distribution  shown  in  Fig.  (6).  The  differ^ 
ence  between  the  patterns  is  not  striking.  The  lb  pattern  is  not  as  intense 
at  the  hot  spot  as  the  la  pattern^  but  t)ilo  was  expected  since  the  lb  pattern 


SECRET 


SECRET 

RESTRICTED  DATA 

ATOMIC  I  NI  lIC.Y  A«  1  1IIS-1 


RM-2115 

29 


had  leoQ  material,  concentrated  near  the  mean. 

A  ccai^TOi'iBon  of  rune  a  and  h  showa  several  Interesting  facts*  The 
largest  cliongc  io  about  30  por  cent  at  the  point  of  maximum  deposition  in 
both  of  the  patterns.  The  log-normal  distribution  of  run  lb  han  fever  par¬ 
ticles  in  the  size  below  than  docs  the  RAND  distribution  of  run  la. 

The  distribution  used  in  run  lb  also  b}iowo  a  higlier  close-in  vfilue  than 
does  the  run  la.  Indicating  too  much  activity  on  larger  particles. 

A  comparison  of  both  of  those  runs  vith  t}ic  fragmentary  data  of  the 
Castle— Bravo  shot  brings  two  facto  to  light.  The  first  point  is  that  both 
the  la  and  lb  distributions  show  fai'  too  large  a  fraction  on  the  islands  of 
the  Bikini  Atoll.  The  second  point  is  that  values  are  low  in  the  northern 
port  of  Rongclap. 

These  errors  may  be  due  to  on  incorrect  wind  analysis,  but  they  could 
elIdo  be  explained  by  a  set  of  incorrect  distribution  functions.  In  any 
event,  the  relatively  slight  change  in  the  patterns  between  an  empirical 
distribution  and  a  mathematical  model  seems  to  be  slight.  This  suggests 
that  an  attempt  could  be  made  to  fit  the  best  possible  log-normal,  or  any 
other  useful  function,  to  the  existing  data. 

Runs  Ic  and  Id  were  made  to  note  the  conporicon  of  the  old  model  vith  the 
new.  Run  Ic  was  the  musliroom  portion,  and  run  Id  was  the  stem;  together 
they  reproduce  a  pattern  vith  essentially  the  some  assunptions  as  the  old 
RAND  model.  The  pattern  for  run  Ic,  the  mushroom  portion,  again  shows  re¬ 
latively  slight  differences  from  runs  la  and  lb.  The  addition  of  the  stem 
material  from  run  Id  makes  little  cliongc  in  the  pattern.  The  greatest  con¬ 
tribution  from  the  stem  is  near  g,round  zero,  which  is  already**  too  high,  and 
near  the  perimeter  of  the  pattern  where  tlie  uncertainty  is  greatest. 
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A  rather  drastic  cliange  in  the  distribution  of  activity  vith  size  vas 
made  In  run  le#  The  amount  of  activity  on  particles  greater  than  95d  radius 
vas  halved,  and  this  exceas  was  added  to  the  particles  between  65u 

in  radiuBo  This  nm  iised  the  spmo  npatial  dlbtrlbutlon  as  runs  la  and  lb. 
The  "hotspot"  for  run  le  was  moved  far  away  i>co8  ground  zero  and  was  reduced 
in  Intensity. 
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Since  the  observations  cm  Cactle  Bravo  were  fragmentaxyj  there  la  no 
useful  purpose  to  be  served  by  continuing  variation  of  parameters  on  the 
set  of  vinds  from  this  shot.  TTie  work  done  has  proven  the  fcaslbU-lty  of 
the  appi*oach  and  Jms  verified  tlie  accuracy  of  the  model  within  the  limits 
oi  the  available  data,  it  is  believed  thi\t  the  next  set  of  data  processed 
can  proceed  on  the  anounqitions  tJiat: 

1#  smoXl  changeo  in  the  dlctributlon  functions  will  produce  only 
negligible  changes  in  the  pattern; 

2*  the  wind  patterns  ore  sulflciently  accurate  to  put  limits  on 
the  poBolble  distribution  functions; 

3.  the  gross  effects  of  the  fallout  from  large-yield  shots  are 
not  affected  by  stem  material. 

It  may  be  Inferred  that  most  of  the  mater ietl  is  lodged  on  particles  in  a 
very  narrow  range  of  sizes  and  that  it  is  concentrated  into  a  very  ruirrov 
range  of  heights  at  the  time  of  stabilization* 

The  next  set  of  rims  will  be  made  on  the  Zunl  shot  of  Operation  Redwing* 
An  attenqpt  will  be  made  to  adjust  the  distribution  functions  so  as  to  pro¬ 
duce  calculated  values  of  the  fraction-down  ^Ich  are  within  the  range  of 
error  meaBurement  at  those  stations  for  which  adequate  measurements  are 
available.  This  then  will  be  uGsunied  to  be  the  optimum  model,  and  the  errors 
of  this  model  will  be  assumed  to  be  a  minimum  for  fallout  calculations. 
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A,  ^P'^TOROLOaICAI.  DATA 

Hun  1  rci’erB  iio  l-nnnn  -riji  1  nut  cilcMltitlC!!^  thnt  v^rc  biiccd  cn  tbc  ’.rindc 
at  the  tliEo  of  the  Cactlc— Bravo  event.  The  vlndc  verc  read  from  the  Dean 
and  Ohmstead  charts  at  8  locations.  Table  A1  gives  these  locations  In  terms 
of  N~S  and  E--W  distances  from  i^ound  zero.  It  vae  necessary  to  Insure  that 
there  were  no  three  points  which  could  bo  simultaneously  collnear  and  closest 
to  the  falling  particle,  because  this  eituiitlon  would  invalidate  the  interpo¬ 
lation  1  unctions  •  tiic  winds  wei‘e  i'uiau  at  two  dlfforont  tlineB  (S.D^  jiuurB 
and  “15 *75  hours  from  H  hour).  For  times  longer  than  8.25  hovira  the  wind 
was  assiflned  to  ho  invariant  with  time.  Table  A2  lists  the  winds  at  the 
eight  points  at  two  different  times. 


Table  A1 

Positions  of  Points,  Relative  to  Ground  Zero, 
Where  Wind  Inl’ormatlon  was  a  Machine  Input 


Point 

Distance  Fast  (n  mi) 

Distance  North  (n  ml) 

A 

-  Ih 

50 

B 

136 

50 

C 

286 

100 

D 

1+36 

50 

E 

“  Ih 

“100 

F 

136 

“100 

G 

286 

“  50 

H 

h36 

“100 
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wind.  Values  for  Run  1 
at  3,25  Hours  after  Sbot  Tljne 
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B.  GEOMl-nBIC  DATA 


Run  la 


Tn  Ic  sdjUS't  "the  ^eom^Tui  Ic  aiuLi  l  fmuiuii  iUnculOius  miu  to 

to  the  yield  of  the  Cnstle-Dravo  event,  it  io  necesBory  to  spread  the  debris 
over  a  lai'ger  region  of  space.  Since  the  lnte/-q'alB  of  the  A  functions  must 
equal  1,  the  actual  activltioc  must  bo  reduced  accordingly.  According  to 
Kellogg^^^^  the  clou/i  diameter  of  n  '5*^  WF  device  is  ~20  n  mi  and  a  15  WT 
device  ia  ^  P9  n  mi*  Since  th^^  Zuni  distribution  of  A^(R)  reduced  to  nearly 
zero  at  ~  17  n  ml  (see  Fig.  8),  it  vns  scaled  up  to  n  ml  for  the  Castle^ 

Bravo  shot.  The  activity  A^  was  therefore  reduced.  In  a  simileir  way  the 
height  of  the  cloud  was  increased,  and  the  activity  fraction  vas  decreased. 

All  of  the  values  of  A  wore  tabulated  to  such  a  degree  that  linear 
interpolation  caused  less  than  1  per  cent  in  the  difference  between  the 
curves  and  the  interpolated  values.  The  tabulated  distribution  functions 
ore  shown  in  Table  A^. 

The  machine  output  gives  the  fraction  of  device  per  square  kllofoot 
at  a  number  of  pre--chosen  points.  Hie  results  of  this  calculation  for 
run  la  are  shown  in  Fig.  Al.  The  isopleths  that  ore  drawn  on  Fig.  9i 
subjective  estimates  of  the  lines  of  equal  fraction  down. 


Run  lb 

The  conditions  for  run  Jh  are  the  same  as  for  run  la  except  for  the 
distribution  of  activity  with  particle  size.  The  distribution  used  in 
run  lb  is  the  log-normal  distribution  with  the  some  mean  and  standard  devlor- 
tion  as  the  distribution  in  run  la.  Table  Ah  shows  the  values  of  as 

a  function  of  r  and  also  the  differences  between  A^(f)  ^or  run  la  and 
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•Xu\}ln  A:^ 

Activity  Function  for  Itun  la 


V.  fVf•^ 

A  /  U  > 

A  /  1  \ 

r 

A  \  1  .TV  4 

1  J■\^W 

1 _ 

yV,  f 

J. 

1  rv V'i  liu.  i 

n 

0 

0.  ^3 

0.15 

0 

.102 

5 

.00024 

0.2^ 

.822 

H 

CD 

10 

.OOO6O 

15.0 

0.33 

1.645 

.174 

15 

.00140 

PO.O 

0.41 

2.46? 

.PJLL 

P.0 

.00240 

0.06 

^.289 

.260 

25 

,00377) 

31=0 

0.93 

4.112 

.316 

30 

.00544 

3O.0 

1.19 

4.934 

.385 

33 

.007<0 

49.0 

1.37 

5.757 

.456 

40 

.01040 

60.0 

1.42 

6.579 

.688 

45 

.01464 

77.5 

l.tl 

7.401 

.764 

50 

.02200 

95.0 

0.68 

8.224 

.495 

55 

.03560 

123.0 

0.43 

9.046 

.372 

60 

.04904 

151.0 

0.27 

9.868 

.303 

65 

.03624 

193.0 

0.19 

10.691 

.252 

70 

.00280 

239.0 

0.20 

11.513 

.218 

75 

.01456 

309.0 

0.16 

12.336 

.191 

80 

.00976 

579.0 

0.11 

13.158 

.169 

05 

.00652 

!  ^»89.0 

0.06 

13.980 

.151 

90 

.00384 

600.0 

0 

14.805 

.131 

95 

.00208 

15.625 

.117 

100 

.00128 

16.447 

.102 

17.270 

.089 

18.092 

.076 

18.914 

.065 

19.937 

.054 

20.559 

.044 

21..  382 

.033 

22.204 

.023 

23.026 

.014 

23.849 

.005 
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for  run  lb.  Figure  A2  bhowa  tlie  I’nllotit  computatlcn  for  run  lb  with  leo- 
pletha  Bubjoctlvely  drfiwn. 

Tabic  Ah 

A^(r)  for  llun  lb  tin  a  Function  of  (r), 
and  the,  Difference  between  lor  Ibon  la  rr.d  for  Hi’Ji  lb 


Index 

r 

A^(r) 

,b 

h.-\ 

1 

9.5 

.2O0 

-.13 

6 

IS. 25 

.389 

-.15 

11 

15.0 

.520 

-.19 

so.o 

.648 

-.24 

21 

24.0 

•  1^1 

-.11 

26 

31.0 

.865 

.07 

31 

33.0 

.918 

.21 

49.0 

.920 

.44 

hi 

60.0 

.883 

.54 

h6 

77.5 

.800 

.61 

51 

rir  r\ 

yj  .V/ 

^Qc 

•  j 

.00 

56 

1^?3.0 

.559 

-•13 

61 

151.0 

.427 

-.16 

66 

195.0 

.310 

-.12 

71 

239.0 

.214 

-.01 

76 

309.0 

.140 

.02 

81 

579.0 

.0O6 

.02 

86 

489.0 

.051 

.01 

91 

600.0 

.028 

.03 

Run  Ic 

In  order  to  be  able  to  check  the  effect  of  the  A(R)  and  A(H)  distribu¬ 
tions  a  model  was  set  which  was  similar  to  the  model  discussed  in 
reference  (1).  The  exponential  decrease  with  height  was  assumed  for  the 
isushroom  and  the  -inv/ivlant  d  i  r.trlbution  with  R  was  retained.  However  97 
per  cent  of  the  activity  was  put  Into  the  mushroom,  instead  of  90  per  cent 
as  In  the  earlier  model,  and  the  remaining'  3  per  cent  was  used  to  make 
run  Id,  The  input  values  are  sliown  in  Table  A5  and  the  final  pattern  is 
shown  in  Fig.  (A3). 
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A, (10  n 

h 

A2(H) 

9.5 
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0 

95 

0.0084 
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0.0067 

Run  Id 

This  run  vao  merely  the  stem  portion  of  run  lc«  The  total  fraction 
placed  in  the  stem  vaa  only  3  cent,  and  the  additioii  of  tills  fallout  to 
run  Ic  made  only  minor  clianf^eo  in  the  pattern.  In  order  to  get  the  gross 
effects  of  the  fallout,  it  ai^pcaj's  unncccscary  to  be  concerned  with  the 
stem.  In  fact  even  this  small  amount  of  stem  fallout  seemed  to  increase  the 
error  in  the  pattern  in  those  firoao  wherein  it  was  detectable. 
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Run  le 

This  confutation  was  doelfTiod  to  nvikr-  n  more  radical  chaii(^e  in  tlio 
dlstrlbxxtlott  function  In  order  to  produce  a  greater  clumr^o  In  the  pattern. 
Since  the  previous  attempts  Indicated  that  t)icre  wao  too  much  dehrla  tno 
close  to  ground  zero,  the  activity  on  the  larger  porticlca  was  drastically 
reduced.  The  activity  size  distribution  of  run  le  is  shown  in  Table  a6  and 
the  resulting  pattern  la  shown  In  Fig.  (A4).  It  may  be  noted  that  the 
activity  Is  spread  further  down-wind  tlian  in  the  otlior  runs,  although  the 
maximum  has  been  considerably  lowered. 


Tvahlc  a6 

The  Function  A^(r)  for  Run  lo 
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Cml  Off. 

Commanding  General,  Elgiith  United  States  Arn^y,  APO  301,  1 

San  Francisco,  California,  ATTN:  ACofS  03 

Commanding  General,  U.S.  Array  Europe,  APO  ^i03,  Nev  York,  N.Y.,  ATTN:  2 

OPOT  Div.,  Combat  Dev  Branch 

Coramandlng  General,  U.S.  Airay  Pacific,  APO  958,  San  Francisco,  Calif.,  2 
ATTN:  Cna  Off 


Commandant,  Command  and  General  Staff  College,  Ft.  Leavenworth,  Ksuisas,  2 
ATTN:  ALrJI5(AS) 

Cccnmandoiit,  Ai*ray  War  College,  Carlisle  Barracks,  Pa.,  ATTN:  Library  1 

Cemmanding  General,  Airay  Medical  Service  School,  Brooke  Army  1 

Medical  Center,  Ft.  Houston,  Texas 
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Director,  Special  Weapons  Develoi)mont  Office,  Headquarters  CONABC,  1 

Ft.  Bliss,  Texas 

Cocp^^dant,  Chemical  Corps  School,  Chemical  Corps  Training  Command,  1 

Ft.  McClellan,  Alnbamn 

Commanding  Crcneral,  U*S.  Am^y'  Chemical  Coi’ps,  Research  and  Devolopiaent  2 
Command,  Washington,  D.C. 

Command  jjig  G,xneral,  The  Fnglneer  Center,  Ft.  Dclvolr,  Va.,  ATTN  I  1 

Asot  Orvlt,  Eng  School 

Commanding  Officer,  Engineer  Recearcli  fuid  Dcvelppmont  I.aboratory,  1 

Ft.  Belvoii’,  Va.,  ATTN:  Ch,  Tecli  lut  Branch 

Commanding  Officer,  Chemical  Wai*farc  Laboratories,  Array  Chemical  2 

Center,  Mfxryland,  ATTN:  Tech  L^br^l^y 

Direetcer,  Technical  Documents  Center,  Evans  Signal  I^ab.,  Belraar,  1 

New  Jersey 

Director,  OperatloriS  Research  Office,  Johns  Hopkins  University,  1 

6935  Arlington  Rd.,  Bothseda  It,  Mm^yliuid 


Commanding  General,  Quartorraastcr  H  and  D  Command,  Quartermaster  R  and  D  2 
Center,  Natick,  Moss.,  ATUi:  CBR  Liaison  Officer 
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NAVY 

ADDRESSEE  NO  OF  CY3 

Chief  of  Naval  Operations,  D/N,  Washington  25,  D.G.,  ATPN:  OP^36  2 

Chief,  Bureau  of  Medicine  and  Surgery,  D/N,  Washington  25,  D.C.,  £ 

ATTN:  Special  Weapons  Defense  Division 

Chief,  Bureau  of  Ships,  d/N,  Washington  25,  D.C.,  ATTN:  Code  346  2 

Chief,  Bureau  of  Yards  and  Docks,  d/N,  Wtuihin^^^ton  25,  D.C.,  ATTN:  1 

Code  D  44o 

Chief,  Bureau  of  Supplies  and  Accounts,  d/n,  Washington  25,  D.C.  1 

Chief,  Bureau  of  Aeronautics,  d/n,  Waohliqi^on  25,  D.C.  2 
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ATOMIC  I  Nl  IJI.Y  A(  I  iOR4 

AFsw)M)if;;rj<TnT./rioN  list 

NAVY  (corrr'D) 


NOjg^CYB 

Chief  of  Naval  Research;  D/N;  WaRhingtou  D.C.,  ATfRi  Ccxle  8ll  1 

President,  U.S.  Naval  War  College,  Newport,  Rhode  Island  1 

Superintendent,  IJ.8.  Naval  Postgraduate  School,  Monterey,  California  1 

Commanding  Officer,  U.S*  Naval  Radiological  Defense  Lab.,  ^ 

Ban  Francisco,  California,  ATTN:  Tech  Info  Divleion  ' 

Commandei^ln-Chlef  pacific.  Pearl  Harbor,  T.H.  1 

AIR  FORCK 

APDRFSSFF  NO  OF  CY8 

Asst  for  Atomic  Energy,  Hq  U3AF,  Washington  25>  U*Cc,  ATTN:  DCS/O  1 

Director  of  Operations,  Hq  USAF,  Washington  25,  D*C*,  ATTN:  1 

Operations  Analysis 

Director  of  Plana,  Hq  USAF-  Washington  25,  D»C«,  ATTN:  Wax*  Plan^  1 

Division 


Director  of  Research  and  Development,  DCS/D,  Hqs  USAF,  Washington  25,  1 

D.C.,  AITN:  Combat  Componento  Dlv 


The  Surgeon  General,  Hq  USAF,  Washington  25,  D*C*,  ATTN:  Bio  Def  1 

Br  Pre  Med  DJv 

CoDunandei^iri-CJiief,  Pacific  Air  Forces,  APO  953,  Son  Fremcisco,  1 

California,  ATTN:  FFCIR-MB,  Base  Recovery 

Commander-ln-Chlef,  Strategic  Air  Command,  Offutt  AFB,  Neb.,  1 

ATTN I  Inspector  General,  Special  Weapons  Branch,  Inspector  Div. 

Commander  Tactical  Air  Command,  Langley  AFB,  VA.,  ATTN?  Doe  1 

Security  Branch 

Research  Directorate,  Hq  Air  Force  Special  Weapons  Center,  2 

Kirtland  AFB,  New  Mexico,  ATTN:  Blast  EFFS  Research 

Commander,  Air  Research  and  Development  Cornmnnd,  Andrews  AFB,  1 

Washington  25,  D.C.,  A'XTN:  RDDN 
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AFSWP  DIf3TJaniJri0N  LIST 
AIR  FOHCK  (C0rrf*D) 


ADDRESSEE  KO  OF  CYS 

Director,  Air  Univereity  Library,  Maxwell  AFB,  Alabaaa  2 

Commandant,  Air  Force  Bcbool  of  Aviation  Medicine,  Randolph  AFB,  2 

Texac 

Consifmder,  /I'  "«••••••  Idge  Reoearch  Center,  L»Q«  Hanscom  Field,  2 

Bedford,  Maas,,  AITN:  CRQ3T“2 

Commander,  AF  Specieil  Weapons  Center,  Kirtland  AFB,  New  Mexico,  3 

ATTN :  Library 


arin-3^  pod  AcrivrriKs 

ADDRESSEE  NO  CF  CYB 

Director,  Weapons  Systems  Evaluation  Groiq),  OSD,  Rm  2E1006,  1 

The  Pentagon,  Washington  25,  D*C» 

Camnander,  Field  Command,  AmiGd  Forces  Special  Wea,pons  Project,  1 

P.O*  Box  5100,  Aj-buquerque,  New  Mexico 

Commander,  Field  Comzaand,  Armed  Forces  Special  Weapons  Project,  1 

P.  0.  Box  5100,  Albuquerque,  New  Mexico,  ATTN:  Tech  TNG  Groiq) 

Chief,  Armed  Forces  Special  Weapons  Project,  Washington  25,  D.C*  15 
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Defense  Threat  Reduction  Agency 

45045  Aviation  Drive 
Dulles,  VA  20166-7517 


CPWC/TRC 


May  6,,  1999 


MEMORANDUM  FOR  DEFENSE  TECHNICAL  INFORMATION  CENTER 
ATTN:  OCQ/MR  WILLIAM  BUSH 


SUBJECT:  DOCUMENT  REVIEW 


The  Defense  Threat  Reduction  Agency's  Security  Office 
has  reviewed  and  declassified  or  assigned  a  new 
distribution  statement: 


-  AFSWP-1069,  AD-341090,  STATEMENT  A 
^DASA-1151,  AD-227900,  STATEMENT  A 
-DASA-1355-1,  STATEMENT  AOf< 


DASA-12  98, 

AD-285252, 

STATEMENT 

A  ^ 

-  DASA-12  90, 

AD-444208, 

STATEMENT 

A-" 

-  DASA-1271, 

AD-276892,*’ 

STATEMENT 

A‘^ 

>  DASA-1279, 

AD-281597, 

STATEMENT 

A^ 

DASA-1237, 

AD-272653, 

STATEMENT 

Al^ 

^  DASA-1246, 

AD-279670, 

STATEMENT 

Ai/^ 

-  DASA-1245, 

AD-419911, 

STATEMENT 

A 

-  DASA-1242, 

AD-279671,‘' 

STATEMENT 

A  ^ 

-  DASA-1256, 

AD-280809, 

STATEMENT 

A^ 

-rr  DASA-1221, 

AD-243886," 

STATEMENT 

"  DASA-1390, 

AD-340311, 

STATEMENT 

-  -DASA-l^Sr 

-AD-717-e97, 

STATEMENT 

A  OK 

--DASA-1285- 

5,  AD-44358-9,  STATEMENT 

' DASA-1714, 

AD-473132, 

STATEMENT 

Ai^ 

-DASA-2214, 

AD-854912, 

STATEMENT 

A*^ 

-DASA-2627, 

AD-514934, 

STATEMENT 

A^ 

^DASA-2651, 

AD-514615, 

STATEMENT 

A 

--EaSA:==-2  536,- 

AD=876697-r- 

^TAtTEMENT 

■A. - 

-DASA-2722T-V3,  AD-518506,  STATEMENT  A 
-DNA-3042F,  AD-525631,  STATEMENT  A 
'DNA-2821Z-1,  AD-522555,  STATEMENT  A' 


